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We calculate the one loop graviton vacuum polarization induced by a massless, nonminimally 
coupled scalar field on Minkowski background. We make use of the Schwinger-Keldysh formalism, 
which allows us to study time dependent phenomena. As an application we compute the leading 
quantum correction to the Newtonian potential of a point particle. The novel aspect of the calcu- 
lation is the use of the Schwinger-Keldysh formalism, within which we calculate the time transients 
induced by switching on the graviton-scalar coupling. 



I. INTRODUCTION 



It has been known for a while that general relativity is, when viewed as a quantum theory, nonrenormalizable. 
Indeed, when gravity is coupled to scalar matter, new local counter-terms are required to cancel the divergences 
^jpj' occurring at one loop ^1], while in the case of pure gravity, new counter-terms occur at two loops [l, That 
does not necessarily mean that quantizing gravity and studying perturbative corrections is meaningless, as long 
as one is interested in low energy effects, which are much below the Planck scale. This effective field theory view 
■ on quantum gravity has been fruitful, and has lead to useful results Observable effects occur primarily in 

cosmology, where as a result of tree level quantization of gravity during inflation, one is led to the Newtonian 
potentials 0, Q that leave observable imprint on the relic cosmic microwave radiation, as well as the seeds for 
the large scale structure of the Universe. The question whether perturbative quantum grav itational effects can 
change the basic tree level predictions of inflationary cosmology is still hotly debated [ol. Il0l| . 
' \ On the other hand, there are well-known tree level quantum effects occurring in curved space-times, the most 
(<~») ■ famous one being the Hawking thermal radiation generated by black holes. There is even a flat space effect 
^-H ' - the well known example being the Unruh effect - which constitutes thermal radiation seen by accelerating 
T-H , observers in Minkowski space. Neither Hawking nor the Unruh effect have so far been confirmed experimentally, 
sjjl ■ The quantum corrections to these tree level effects are as yet not well understood. In particular, it would be of 
• ^ , interest to investigate whether such quantum corrections can induce a significant back-reaction on the background 
' space-time that would correct some undesirable features that occur when fields are quantized on black hole spaces, 
H , such as quasi-normal modes (which essentially tell us that it is inconsistent to quantize massless fields on black 
■ - - ■ hole backgrounds). The main message from these studies is that: we shall be able to incorporate self-consistently 
the back-reaction from fluctuating quantum fields on black hole backgrounds only if we allow space-times to 
become dynamical. In other words, a self-consistent semiclassical gravitational theory will require a (quantum) 
modification of the Birkhoff theorem (whose classical version states that a spherically symmetric distribution of 
matter can induce static, radially dependent, metrics only). 

In this paper we do not address this very interesting question, but we point out that the Schwinger-Keldysh 
formalism |ll| - [l5| is suitable for such studies. As an example of how to apply the formalism, we provide a 
perturbative calculation of the quantum correction to the Newtonian potential generated by a massless scalar 
field. This represents a baby version of the more interesting problem of quantum corrections to black hole 
space-times. Our treatment generalizes the work of Park and Woodard T^, in that we consider not just the 
minimally coupled massless scalar field, but also include a nonminimal coupling of the scalar field to the Ricci 
curvature scalar. Another important difference is in that Park and Woodard assumed that gravity was turned 
on at ^0 ~^ —oo, thus preventing any time transients from occurring. We, on the other hand, assume that the 
scalar field decouples from gravity at early times (t < to), and that the coupling turns on at t = Iq, where to is 
some finite time. In practice, this can be realized, for example, by a Higgs mechanism, which gives a large mass 
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to the scalar at t < to, and the mass vanishes at t > fp- Such a scalar will effectively decouple from gravity 
at t < to and on the length scales L > l/m^, and will induce a nontrivial gravitational effect for t > to- In this 
paper we demonstrate that the Schwinger-Keldysh formalism allows us to follow in a manifestly causal manner 
the gravitational field induced by such quantum scalar field fiuctuations. 

Of course, our work is not the first to address quantum corrections to Newtonian potentials. Such studies were 
pioneered by Donoghue [13, [13, and followed by many others [l9| , l23l . with differing results. The problem was 
again revived in the 2000s by Bjerrum-Bohr et al. [13) Butt [24| . Faller [2^ etc. where also some of the 
graviton vertex corrections were considered. Also the works of Dalvit, Mazzitelli, Satz and Alvarez [2^[23l are of 
interest to us. Indeed, in Appendix C we show that there is a rather intricate connection between the methods 
used in Refs. [H, [23| and our methods. 

The paper is organized as follows: in Sec. II we briefiy show how to obtain the classical Newtonian potential. 
Section III is devoted to calculation of the one-loop graviton vacuum polarization induced by a nonminimally 
coupled scalar field, while Sec. IV presents an application of the main result of Sec. Ill: the quantum one-loop 
correction to the Newtonian potential. In Sec. V we summarize our main results. In the Appendixes we review 
the Schwinger-Keldysh formalism (Appendix A); some subtleties in deriving the Newtonian potential of a point 
particle (Appendix B); we show how to expand the 2PI scalar bubble diagram to get the stress-energy tensor and 
the graviton vacuum polarization tensor (Appendix C); and we present details of the one- loop vacuum polarization 
calculation (Appendix D). 

Unless stated explicitly, we work in natural units where h = 1 = c. 



II. CLASSICAL NEWTONIAN POTENTIAL 



We are interested in the gravitational response of a static, point-like particle of a mass A'l in the particle's rest 
frame, with the classical stress energy tensor in Minkowski background 

T(^;=AM°<5«5^-i(f), (1) 

where D denotes the dimension of space-time. For gravity we take the classical Einstein-Hilbert action: 

Seh = J d^xy/^R (n^^lGnGN), (2) 

where R is the Ricci scalar, g is the determinant of the metric tensor g^i,, and Gn is the Newton constant. 

In order to get the classical Newtonian potential, we first expand the metric tensor around the fiat Minkowski 

metric 77^1,: 

gt,u ^ ^l^iu + h^i, , ry^i, = diag(-l, 1^^), (3) 

D-l 

where, in general, h^u — h^v{x), x = (t,x). The full classical gravitational response to a point mass is obtained 
by varying the action ^ and setting it equal to the classical stress-energy tensor, SSm/Sh^^" = — -y/^/2T^i.. 
Here we are primarily interested in the leading gravitational response (Newtonian potential) to a point particle 
at rest, whose stress-energy tensor is given by ([1]). To accomplish this we need the Einstein-Hilbert action to 
quadratic order in /i^j, (for detailed calculations see Appendix B): 



Seh = TT^ I d^x 



h^'^'L^^p^hP" + 0{hl^) , (4) 



where L^^^pa- stands for the Lichnerowicz operator in Minkowski background: 

Lpupa = d(pri„)(pd^) - ^rjpi^prj^^^d'^ ~ \^{'r]yivdpd„ + T]pcrdpdt,) + ^7jf,uVp<y9'^ ■ (5) 
For a later use we rewrite the Lichnerowicz operator as a combination of two simpler operators: 

LpupcF — Dpup^ "^"^pLv-Dpa , (6) 
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where 

Dpa = dpda - ?7p<t9^ , Dp^pa = d{pr]„){p,dy) ~ ^Vt^(pVa)iyd'^ - ^'npadf.d^ . (7) 

The Lichnerowicz operator ([5]) possesses the following symmetries: it is symmetric under the exchange of the first 
two indices (/i ^ i^); of the last two indices (p O a), as well as under the exchange of the first two and the last 
two indices {fiv O pa). 

Now varying the action Q with the stress energy tensor ([T]) yields 

y '5°5°M5^-i(f). (8) 

The same equation can be obtained by linearizing the Einstein equation around the Minkowski background, which 
is a simpler procedure. As we show in Appendix B, the Newtonian potential in the longitudinal (Newton) gauge 
Eq. ^ gives the solution: 

4o^(x) = 4"^(x) = ^, (^ = 1,2,3), (9) 

where the superscript (0) emphasizes that the metric components refer to the classical solution. In the 
following two sections we show how to use this solution to construct the quantum corrected one-loop Newtonian 
potential. 



III. GRAVITON VACUUM POLARIZATION 



Here we derive the leading quantum contribution to the graviton vacuum polarization tensor due to a nonmin- 
imally coupled massless scalar field if, with the action: 

S^^J d^x^S^-^{dp^){d,ip)g'^'' - ^^R^'^ , (10) 

where ^ measures the coupling strength of to gravity through the Ricci scalar R. In this paper we focus on 
the quantum corrections around Minkowski space, hence we expand the metric tensor around the flat Minkowski 
metric (c/. Eq. ©): 

9 pi, = Vf^iy + i^hp^ . (11) 

Since we are here primarily interested in the one-loop graviton vacuum-polarization, in (jlip we have extracted 
the gravitational coupling constant k defined in Eq. 1^, as can be used as the loop counting parameter of 
perturbative gravity. Note that hp^{x) in Eq. ([TT|) can be understood to contain both classical and quantum 
contributions to the Newtonian potential, which is in the spirit of the rest of the paper. One can write Eq. pip 
in the more general form, g^i, — g'^p} + nhp^, where g^i) denotes a classical background metric around which 
one expands. In the context of this work, the natural choices for g'^p}, are the metric of the classical Newtonian 

potential, g^p} = ripi, + h^p} with the elements of given in Eq. ([9|), or by the Schwarzschild black hole metric. 
Since choosing these background metrics would entail significant technical complications, we leave their treatment 
for a future work. 

In our model ^ and (jlOp the graviton acquires quantum contributions both from the graviton and scalar quantum 
fluctuations. To calculate the graviton one loop vacuum polarization one needs the cubic and quartic graviton 
vertices, which can be quite straightforwardly extracted from Eq. (|82p . For simplicity in this paper we focus on 
the one- loop scalar contribution to the graviton self-energy (vacuum polarization). To this end we need the cubic 
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and quartic scalar-graviton vertices, which can be easily extracted from the scalar action (jlOp : 



(12) 



d^x \ h'^'ix) 



X {daip{x)){dpip{x)) 

h 



(13) 



where we made use of Eqs. ([MI), (El) and ([SHZl). 

In the spirit of the Schwinger-Keldysh formalism, cubic and quartic parts contribute to the interaction action 
as: 



SintVp^ ,h+p,ip ,h^p]= Sh^^ltp^ ,h+p]~ Sh^vVp + Sh2^2[Lp+ ,h+p]- Sh2^2[ip ,h^p\ 



(14) 



These interactions can be used to generate the graviton tadpole and the graviton vacuum polarization induced by 
the one-loop scalar fluctuations. The scalar diagram that contributes to the graviton tadpole is shown in Fig.m 



Vh 



FIG. 1: The graviton one-loop tadpole. Only the scalar loop contribution is shown. The scalar propagator A is the blue 
dashed line, while the graviton insertion h^v is the red solid wavy line, and k = \/167rG]v is the gravitational coupling 
constant. 

The three scalar field diagrams that contribute to the one-particle irreducible (IPI) graviton vacuum polar- 
ization at the one-loop order are shown in Fig. [2l The first diagram is the local contribution generated by the 




+ 



K 



5Z 



FIG. 2: The scalar diagrams that contribute to the graviton one-loop vacuum polarization. 5Z stands for the one-loop 
counter-terms. 



quartic action ([T^. the second is the non-local contribution generated by the cubic action (IT^ squared, and the 
third is the counter-term. 

In this paper we work with the IPI effective action, which consists of the tree-level (Einstein-Hilbert) part, the 
matter free action, and the one-loop contributions with zero, one, two, three, etc. graviton insertions. The one 
loop bubble diagrams with no graviton insertions are shown in Fig. |31 Formally, we can write the one-loop IPI 
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FIG. 3: The scalar and graviton one-loop bubble diagrams. puA^,^ denotes the graviton propagator. 



contributions to the effective action as, 

ri[V,A] = -^Trln[A''n2:;x)]-^Trln[^.A»»(x;x)]-i ^ f d'' x^'^ {x)a6'^%])'^''{x) 

a,b=± •' 

- I E h""^! d'^x'K^ix) [/,n^J {x;x')K\x') (15) 

„ 1, l_ ^ 



OO ^ 

E IT 



a,fc=± ' 

^ / d^x.Kl-^^ {x^) [ d^'x^Kr (^2) • • •/ d'^x,,Kf- (a^n)V^;,V?p.:?..^„.„ a:^; . . . ; 



n— 3 ai,a2,..,a 

where A°''(x; x') denotes the scalar Schwinger-Keldysh propagator, which for a massless scalar field on Minkowski 
space can be found in Eqs. (|5Dl - f5^ in Appendix A, Tr denotes both a trace over the space-time indices and an 
integration over space-time, tI^J°'°'{x) denotes the tadpole in Fig.[Tl [ ^^11^^] (x; x') is the graviton vacuum polar- 
ization tensor, for which the contributing one-loop diagrams are shown in Fig. [21 and y^^viii2V2 ^i^vn (^i' -^2; ■ • • ; Xn) 
stand for the one-loop vertex corrections with n > 3 graviton insertions. An example of a vertex correction dia- 
gram is shown in Fig. 21 In order to get the complete one-loop correction to the graviton, one needs to evaluate 



FIG. 4: A one-loop vertex diagram with one scalar loop and six graviton insertions. 

all of the diagrams contributing to the one-loop effective action p^ . This is hard. However, the corrections with 
more and more graviton insertions are expected to be progressively smaller and smaller, and hence it is often 
enough to calculate the leading non- vanishing correction. In this paper we evaluate this leading order correction, 
which corresponds to two-graviton insertions. In fact, all leading order quantum effects on the graviton are 
embodied in the graviton vacuum polarization tensor, and the contributing diagrams are shown in Fig. [51 

The effective action (|T5)) can be derived by a Legendre transform of the free action. Here we sketch a heuristic 
derivation of (jlSp . which misses the scalar and graviton bubble diagrams. The loop contributions to the effective 
action can loosely be defined as an expectation value of e*'^'"', where S-mt is the interaction action given in (fTl]) : 

e^r, ^ ^g.s.„.^ ^ 1 ^ ^^^+ ^ _ s-^^^ ^ ^^^+^^ _ s-^^^^ _ 1 ^(^+^^ _ ^-^^)2^ ^ 0({h%f) , (16) 

where (•) denotes an expectation value with respect to a chosen state which in our case is the Minkowski 
vacuum, corrected for the fact that the interaction with the graviton switches on at a finite time to. This 
difference allows us to investigate the thereby induced time transients, and in particular how causality affects 
and limits the growth in time of the terms induced by quantum effects. In Eq. (|TB| we introduced the notation, 
^hifif ~ ^h.ipipVP^ , h^v\ ^iid <S'^2^2 = SY{2^p2[^p^ ^ h^v]- ^his work we neglect the vertex corrections to Fi, which 
occur at the order 0{h^^) and higher. It would be of interest to investigate whether such corrections change any 
of the results presented in this work. For a discussion of the role of vertex corrections to the Newtonian potential 
see Refs. [U The first non-trivial term on the right-hand-side of (|T6| yields the one-loop graviton tadpole 
of Fig. [1] and thus the scalar one-loop contribution to the stress-energy. The second contribution yields the local 
one-loop contribution to the graviton self-energy, shown in the first diagram of Fig. [2l The disconnected part of 
the third contribution is the tadpole squared, which is a part of the geometric series that is needed to reconstruct 
exp(iri); while the connected part of the third contribution represents the nonlocal contribution to the one-loop 
self-energy, diagrammatically shown by the middle diagram of Fig. [21 etc. 
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More economical is the two-particle irreducible (2PI) formalism, in which one formally writes the effective 

action in terms of scalar propagators A(a;; cc'; [g^^J]) in a general background g^S- The one-loop contributions to 
the 2PI effective action are then the scalar and graviton bubble diagrams from Fig. |3l 



ri%[A, ,.A,„] = -^Trln[A-(2/;y; [g^])] „ ^Trln[ ,,A;:iy;y; g^'})] , (17) 




FIG. 5: The two-loop diagrams contributing to the 2PI effective action. 



the first term in Eq. ([TTl) with respect to gfiv , and insert this into the classical linearized equation for the graviton 
with a static point mass source to study the leading quantum correction to the Newtonian potential, which in 
that approach originates from the quadratic term in the expansion of the scalar propagator A"'' in Eq. (jl7p in 

powers of ft.^^ = — t^^j/. In Appendix C we show that that approach must lead to the same answer for the 
quantum corrected Newtonian potential as obtained from the one-loop graviton vacuum polarization, which is 
the approach advocated in this paper. A detailed comparison shows that the results indeed agree, representing a 
non-trivial check of our work. 

Let us now go back and discuss various contributions to Eq. (jlSp . When the scalar and graviton propagators 
are calculated on Minkowski vacuum, they do not depend on and hence the one- loop bubble terms in (|15|) 
contribute to the effective action as a constant, and can be neglected. Therefore, the one-loop graviton effective 
action in our theory has the following form (c/. also Eq. (|113l) in Appendix C), 

Ti[h^.] = SEH[ht.]-l J d^^/ir(^)TW±(x)-i J d^xd^x'h^^^ix)[^^U%]{x;x')h^^^{x') + 0{ih^,f). (18) 

We shall first argue that the tadpole and the local contribution to the graviton vacuum polarization vanish and 
then evaluate the nonlocal contribution to the one-loop graviton vacuum polarization. 

A. Local contributions 

Since we work in dimensional regularization, the local diagram ~ which yields a quadratic divergence - is 
automatically subtracted, and hence will not contribute. To see this in some detail, from Eq. ((T6)) we infer that 
the tadpole contribution to the effective action yields 

Ttadpolo = ^h<p<p[h±'] — {Shipip — S,^^^) , (19) 

where S^^^ are defined in Eqs. (IT^ and p^ . If we take look at the structure of S^^^ in Eq. we see that 
the contributing terms can be written as 

{ip{xf)=iA++{x;x) and {da(p{x)di3ip{x)) = dad'piA++{x]x')\x^x' ■ (20) 

The first term in Eq. ([20)1 equals the coincident {x = x') Feynman propagator «A++ cx l/Aa;;^^^ (see Eq. (1501) ). 
which is in dimensional regularization defined as the analytic extension of the propagator where it is defined, i.e. 
in '^[D] < 2, where iA^^{x; x) = 0. Hence, this term vanishes in the dimensional regularization. In order to 
show that the contribution of the second term vanishes, first note that 9q9^iA++(x; x') — S'^Spi6^{x — x'). This 
implies that when either a ^ or /3 7^ 0, this term vanishes, since in dimensional regularization it equals the 
analytic extension of 9Q,9^«A++(a;; x')\x^x' from 5R[_D] < (where it vanishes) to the whole complex plane. When 
both a = = /3, then this term yields i5^{ X — j-x — which, cru be shown to Vcinish in dimensioncil 

regularization. We have thus shown that the tadpole vanishes in dimensional regularization: 



Ttadpolc — 



(21) 
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Let us next consider the local contribution to the one-loop graviton vacuum polarization: 



r 



local 



s 



(22) 



From the structure of the quartic action we see that the same terms as in Eq. ([20]) contribute to Fiocai, and 
therefore, in dimensional regularization, 



r 



local 



= 0. 



(23) 



Non-local contributions 



As previously elaborated, the only term in Eq. (fT6)) that contributes to the one-loop effective action (|T8)) is 
contained in the Wick-contracted part of the expression iSThipip = — SJ^^^Y): 



6Tn^^[^^,h%] ^ -\j d^x I d^x'h^^ix) [%U%] {x;x')hi^ix') 



(24) 



_ ( ^ ) ( ^ ) / d'^xd'^x'^^"'' 



2 \~) It 

1 

2 



X 



X 



/i5f(x') 



where the primed partial derivatives stand for derivatives with respect to x' . Next, we perform a detailed 
calculation of the vacuum polarization tensor in Appendix D and note that the derivative operators appearing in 
the vacuum polarization tensor (jl28l) can be written in the useful form: 



^ apiiv ^ per ^ iivcxp y ^ pa ^ '/ ^ pa 



^ pa 



(25) 
(26) 



where L^j^^ap is the Lichnerowicz operator ([SHZ]). Making use of the results from Appendix D, we can split the 
graviton polarization tensor (jl28p into a non-local finite part and a local divergent part as follows: 



^ (x; x') = ^ (x; x')+^ [±,n± ] '^'^^ (x; x') 

where the renormalized graviton one-loop vacuum polarization tensor reads 

z[±.n±]('-™\x;x') = -(±)(±): 



(27) 



K 



307207r4 

- n"^ n 

pa '/ ^ pa 



[ln2(^2^x^^) - 21n(^2^x^^)] . (28) 



This is the generalization of the graviton vacuum polarization induced by a minimally coupled scalar (^ — 0) 
calculated by Park and Woodard in jl^ to the case of a massless, nonminimally coupled, scalar (£^ 0), and it is 
one of the main results of this work. When ^ = 0, our result (I28p is in perfect agreement with [l6|, representing 
a nontrivial check of our work. 

The divergent part of the vacuum polarization tensor (|119p can be extracted by inserting the divergent part of 
Eq. into ([^51) . The resuh is 



D\ ,,D-4 



647r^/2(D2 - 1){D -2){D- ■i){D - 4) 



(29) 



X UD, 



D 



D - 1 



8(^^-1) U 



D-2 



As we show in the next section, this term has the right structure such that it can be subtracted by local counter- 
terms, which ought to be added to the effective action to complete the one- loop renormalization program. 
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Renormalization 



In order to cancel the divergent part of the graviton vacuum polarization ((29|) . from the form of the effective 
action (fT8|) we see that it should comprise the counter-term(s) of the form 



(30) 



In spite of the appearance, this counter-term pOI) is local (because of the delta function in Eq. (1^^ ). as it should 
be. 

As we shall now show, two counter-terms that arc quadratic in the curvature tensors are needed to subtract 
the divergences in (|29p . To see this, first note that the linearized form of Ricci scalar and Ricci tensor are (c/. 
Eq. (HI)) 



(31) 
(32) 

where we expanded the metric tensor around Minkowski space, g^^i, = r/^^ + nh^^^. To quadratic order in h^^^ the 
following two counter-terms^ 



<5re.t. = ai 5](±) j d"^gxR[9%f + a2Y.^±) J d''x^R^4g%]R^'^[g%] 



(33) 



and can be recast as 



<5rc.t. = ^(a3)±±«:2 f d°xd''x'h^^^ix)\(a,D^,D'^, + a2D^p^,D'''^^„)5''ix - x')W{x') 
±± •' ^ ^ 

Comparing this with Eqs. (I^^l and (P^ . we can read off the coefficients ai and a2 in p4l) : 



(34) 



ai 



r(f)M^-^ 



1287r^/2(L>2 - 1)(L» - 2){D - 3){D - 4) 
327r^/2(i:>2 - 1){D -2){D- 3){D - 4) 



D 



D - 1 



D-2 
4(^-1) 



(35) 
(36) 



where aj'^"^ and a'^'^^ represent the finite parts of the counter-terms, which are to be fixed by measurements. One 
often uses the minimal subtraction scheme, in which one expands the coefficients ai and a2 around 1/{D — 4). 
In this case their divergent parts are 



Jdiv) _ 



1 



9607r2 



1 



-3+n^-6 



D-4 



and 



(dr 



1 



D-4 



960772 p _ 4) 



(37) 



and - when f = - agrees with Eqs. (3.3) and (3.34) of Ref. Q. ^ This completes our discussion of renormalization. 



D. Retarded Self Energy 

By varying the action with respect to the fields h'^ or h'^ and setting h'^ — — /i^" 

= 0, 



(38) 



^ In fact, Veltman and 't Hooft get a factor two larger result for oi and 02, but they work with a complex scalar field, which can be 
decomposed into two real scalar fields, thus explaining the difference. 
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the quantum-corrected equation of motion for the metric perturbation /i^i^ is 

L^,p.hP%x) + j d^x' [,,,n-*] {x- x')h'"'{x') + 0{{h^,f) = '^SlSlMS'ix) , (39) 

where 

U^Hp.]™' {x; x') = [+,n+ ] (-^^"^ (x, x') + [+,n;,] ^-^^"^ (x; x') (40) 

is the retarded graviton vacuum polarization, which we use in the next section to compute the quantum-corrected 
Newtonian potential. 



IV. THE QUANTUM CORRECTED NEWTONIAN POTENTIAL 



Rather than attempting to consistently solve Eq. (p9)) . we shall solve it perturbatively, i.e. we make the 
following perturbative Ansatz for the graviton field: 

K,{x) = h^^}{x)+h^^){x), (4f) 

where h^^} stands for the classical solution © obtained by solving the classical part © of the full equation of 
motion ([Mt . 

To find the quantum correction we need to solve the perturbative equation 



L^.p.hP^^^ (x) + j d^x' l^.U.^Y''' (x; x')hP'^(''^ {x') = , (42) 

where [^,ynpo-]'°' is given in Eqs. (pSj) and (HU]) . ^ 

In order to evaluate [p.j^IIpo-]'^''', it is useful to split the logarithm function in (|28p into its real and imaginary 
parts: 

ln(/i2Aa;^+) = In + Ar^ + ^e)] = In + Ar^)] + n:Q{M'^ ~ Ar^) (43) 

\n{li^l\x\_) = In + Ar^ - le sign(At))] = In + Ar^)] - n:Q{M'^ - Ar^) sign(At) , 

with the abbreviations Ar = ||a^— x'\\ and Ai = t — t'. The retarded graviton vacuum polarization is then 

Un-](x;x') = -_g-jxL,.„,|i?"V-^/"^^P.[J + 

X e(At2 - Ar2)e(Ai) [In |^2(Ar2 - At'^)\ - l] . (44) 

Recalling that h^°'>P''{x') = 2GnM5P'' /r' (r' = ||f' ||) for p ct = 0, ...,3 and inserting it in Eq. ^ together 
with (|13]) we get 



(45) 



where 



F"'^ = F(r, = / dt' f d^x'e{At - Ar) [in fi^{At^ - Ar^) - l] 1 . (46) 



^ Had we attempted to find an exact solution to II39I I. we would need to solve an integral equation for /i'^" = {2G m M /r)Ho{Kr) 
/i"*^ = 5^^ {2G iq M /r)H\{K,r), where Ho(z) and H\{z) are the sought-for functions, whose asymptotic series (around r = oo) begins 
as H\^2{z) = 1 -|- 0(1/2^). It would be of interest to find out what are the coefficients of this asymptotic series and whether it 
can be resummed, i.e. whether it can be analytically extended to the whole complex z-plane. Understanding this function is 
important, because it would tell us whether quantum corrections can resolve the r = singularity of Newtonian gravity. Pursuing 
this analysis is, however, beyond the scope of this work. 
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Here we are not interested in homogeneous solutions of the operator Lf^^p^r, and thus Lpi,p^ can be dropped out 
from both sides of Eq. (|45|) . by which (|45|) simplifies considerably. 

It is much easier to carry out the integration in ()45|) if we change the variable of integration x' to Ar = x — x' , 
i.e., d^x' = d^Ar = 27r(Ar)^d(Ar) sm6d6. The integration over 9 yields 



sin 9d9 



1 



c?(cos 9) 



r + Ar — \r — Ar| 
rAr 



(47) 



where r = ||x || and we choose r' = zr' . We thus have 
ft 



-1 ^7-2 + Ar2 - 2rArcos(6') 
We thus have 

F{r,t)=ATT dt' d(Ar)(Ar)2{ln[^2^At2-Ar2)]-l}|^e(Ar-r) + ie(r-Ar)|e(At-Ar). (48) 
For simplicity, we decompose this integral into two parts, namely I a and Ib'- 

F = Ia+Ib, Ia = in dt' d{Ar){Ar){\n[^^{At^ - Ar^)]-l}e{Ar-r)e{At- Ar) (49) 

J t o Jo 

/ dt' d{Ar){Ar)^{\n[n^{At^-~Ar^)]^l}e{r-Ar)e{At^Ar). (50) 

J to Jo 



An 



If we now take Ato = t — to and carefully treat the 0-functions, for Ia we get the following result: 



Ia ^ An d{At) / d(Ar)Ar {\n[fi^{Af - Ar^)] - l} e{At - r) 

Jo Jr 

- ^Atl + aOAtor^ - 22r3 - 6r=^ In | ^^^-^ ) + 

V Ato - r / 

- iAto{Atl - 3r2) In [^^^(Aig - r^)] le(Ato - r) 



(51) 



27r 



In order to calculate Ib , we first write it as a sum of two 0-functions 

^At, 



12r^ ln(2/ir) 



47r Z"^*" r 

Ib^— dAt \ d{Ar)Ar^\n[fi^{Af - Ar^)]~l}eiAt-r) 
f Jo VJo 

+ ^ d(Ar)Ar2 {ln[^2(^^2 _ ^^2)j _ _ ^^^j 



(52) 



Note that, after integrating over Ar, the second integral in Ib gets divided in two parts, thus explaining the 
origin of the time transients. Schematically this means 



^ °dAt| ••• |e(r-At) — > °dAt| ••• |e(r-Aio) + j|| dAi| ••• |e(Aio 



(53) 



It is now convenient to break the overall result for F into a sum over two parts, each proportional to one 
G-function: 



F{r,t) = Fi{r,t)Q{Ato-r) + F2{r,t)Q{r-Ato) , 



where Fi and F2 evaluate to 

3 

F^[T,t) = 47r^ 
6 



\ N 25' 
ln(2Mr) - - 



+ 47r- 



Atl - 



6 



(Ato + r)^ 
2r 



ln[/x(Ato + r)] 



(Ato - r) 
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2^ ln[/x(Ato - r)] - y Ato 



F2(r,t) ^ 47r^ 
or 



25 

ln(2/iAto) - — 



One can easily verify that d^F satisfies 

d^F{r,t) = [d^F^{r,t)\ e(Ato - r) + [d^F2(r,t)\ e(r- Ato) : 



(54) 

(55) 
(56) 
(57) 

(58) 
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and analogously for d*F. This simple result holds true because, not only Fi and F2 are continuous at the causal 
boundary Ato = r, but also the pairs {{do + dr)Fi, ido + dr)F2}, {d^Fi,d^F2} and {{do + dr)d'^Fi, {do + dr)d'^ F2} 
are all continuous at Atp = r. Finally, from (|55H57| and (f58| we obtain: 

(59) 



d^F . 16 Ji^e(A.o - .) + i!li^e(. - Ato)l 

r r J 



By inserting this into (|45p . we get the desired result for the quantum one- loop correction to the Newtonian 
potential: 



= ^ {5^ (1 + 2«) - - ^ (-1 + ^ - A,„) + 1^(1 - inr - Aw} (60) 

k^GnM r 1 



where ^ denotes 

|=i+3ofe-^) • (63) 



6 V 6, 

Eqs. (|60H63p constitute our second main result. Apart from the static results, which are reproduced in the limit 
when Ato 00, Eqs. ([60l[62|) contain time transients both on the light cone (where Ato = '') and outside the 
lightcone (where Ato < f). The coincident time divergences oc l/(Ato)^ and oc l/(Ato) in (|5DH5^ outside the 
light cone signal initial time divergences, which occur as a consequence of the sudden switching of the graviton- 
scalar coupling, and can be removed by either resorting to an adiabatic switching ^2§\ , or by a suitably modifying 
the initial state at t = to [sol [30i]. However, performing this initial surface renormalization is beyond the scope 
of this work. The appearance of the time transients in ()60H62p is the main novelty of our approach to quantum 
gravitational corrections on classical space times, and can be of help to understanding the dynamical quantum 
gravitational back-reaction. 

In the limit when Ato 00 (to — ?► ~oo) the second and the third part of and h^^^-^'' in ([501 - 1511) vanish 

and all of /i*''^^) in (15^ vanishes, and one recovers the static result. This is to be expected, since in the limit 
when to — > —00 time transients should disappear. In this limit our results reduce to the ones obtained by the 
in-out formalism, in which calculations are usually performed in momentum space, and no time transients are 
allowed. Indeed Refs. [2^ [l^l have evaluated the one-loop correction for a nonminimally coupled scalar in the 
same (longitudinal or Newtonian) gauge, albeit by different techniques (see Appendix C). Their results are shown 
in Eqs. (31) and (32-33). When compared with ours (I5D1 - I5T|) (with Ato cjo), the coefficients of h'^^'^^^ and 



agree in magnitude, but they have the opposite sign. Comparing, however, with Park and Woodard 16|. 
which calculate the effect of a minimally coupled scalar only, our results agree perfectly. To show that, observe 
first that in the longitudinal gauge we use, and in the limit Ato —5^ 00 (c/. Eqs. (IBOMTI) and ([M])) we get the 
quantum-corrected Bardeen potentials 

^^_G^_ G^l + 2|^^/^x GnM GlMl-21 ^ ^/ 1 

r GOtt r'^ J ^ r GOtt r-^ 



whose leading quantum parts in the ^ = (^ = 1) case are, 

(,) 5^0 _G%M ^0 G%M 

On the other hand, relation (|106p and Eqs. (55-56) of Ref. [16] tell us that, in the gauge used in Ref. [l6|, 
/loopw = G%M/{10nr^), h^pl^ = ~G%M/{10TTr^) and h^plj = -G%M/ {SOnr). Inserting these into Eqs. ^ 
and (fTOT)) we get = -GjfM/{20TTr^) and ^-^^^ = G%M/{60TTr^), which perfectly agree with ([eS]). 

An interesting question is: in which physical situations one expects to see time transients of the kind exhibited 
in Eqs. (|60l[62l) . One such situation would be a very massive scalar field which at time t = to suddenly becomes 
massless {e.g. via a Higgs-like mechanism). Such a scalar would effectively decouple from gravity at t < to, 
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and would begin acting gravitationally at t — Iq. Equations (|60H62p describe a good approximation to the 
gravitational field generated by such a scalar. If the Higgs mechanism would induce a significant change in the 
vacuum energy, one would need to consider the change in the Newtonian potential from a de Sitter background 
to a Minkowski background. While certainly interesting, this calculation would be technically demanding. For a 
recent calculation of the one- loop vacuum polarization in de Sitter background see Ref. [s^- One may object that 
the work presented here is mainly of academic interest. While this is to an extent true, we emphasize that there 
are many similar situations that are realized in nature, where time transients are of crucial importance. One such 
example is a shell of collapsing matter, that eventually forms a black hole. While understanding the quantum 
back-reaction in such a dynamical setting is beyond the scope of the (momentum space) in- out formalism, as the 
present work indicates, it is well within the scope of an in-in treatment. 



V. CONCLUSIONS 



We have calculated the one-loop graviton vacuum polarization induced by a massless non-niinimally coupled 
scalar field by making use of the Schwinger-Keldysh formalism. We have then applied our result to calculate 
the perturbative quantum correction to the Newtonian potential of a point particle. The novelty of our approach 
are the time transients (|60H63p . which naturally appear within the Schwinger-Keldysh formalism. Such transients 
would occur, for example, in the case when the scalar mass would change from a large value to zero by a Higgs-like 
mechanism. 

There are many directions in which the calculation presented in this paper might be developed. One such 
direction is to consider quantum effects due to fermionic, vector and tensor fields. Other possible avenues 
comprise a study of time transients induced by quantum effects in the process of formation of compact stars, such 
as neutron stars, black holes, boson stars or gravastars [s^. In particular, these types of studies would improve 
our understanding of the question of quantum backreaction. For example, growing time transients would indicate 
that the background space-time ought to be modified. 



Appendix A: Schwinger-Keldysh formalism 

In this section we give a short overview of the Schwinger-Keldysh or in-in formalism [Tll - fisj by making a 
comparison with the more conventional in-out formulation of quantum field theory (QFT). Here we shall follow 
mainly Refs. [l3,[il]- Within the in-out method the vacuum-to-vacuum amplitudes {VL, tout|f^, iin)j are calculated 
in the presence of an external source J in virtue of the path-integral representation of the in-out generating 
functional Z[J]: 



(f^,iout|f^,iin).7 = ^[J] = j Pv^exp i{s[ip\+ j d^xJip) 



(66) 



where |r2,tin/out) represent the in and out states in the remote past and future, respectively, S[ip\ is the free 
(classical) action for a (scalar) field '^{x), J{x) is a source current, and Vip = Y\^d'f[x) is a path integral 
integration measure. In flat space-time the in and out states are mostly expressed as a superposition of free plane 
waves defined globally; in curved spaces however the in and out states can differ (in that they are defined only on 
part of the space) and the out states may even not be known. For these type of situations the in-in formalism is 
more suitable. The partition function Z[J] of the in-out formalism generates time-ordered correlation functions 
between different states: 

(5" 

(r!,tout|T{^(a;i)---^(x„)}|f],U = H)"T7r^ — --Z[J]|,;^o • (67) 

In nontrivial gravitational backgrounds the final state of a system usually cannot be predefined, and hence one 
resorts to calculating expectation values with respect to an in state. In doing so one first defines the in-in 
generating functional by means of the two external sources J_|_ and J_ : 

Z[J_, J+] =j_ {il,tin\^,tin)j_^_ ^^{^, tin\a, tout) J^{a,tout\^, tin) J+, (68) 



where the sum goes over a complete set of out states. While in the standard in-out formalism one evolves the in 
state in the presence of a source J and compares it with the out state, in the in-in framework one evolves the in 
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state in the presence of the two different sources J+ and J_ and compares the resuhs in the remote future. In 
other words, expression (j68p can be understood as the m-vacuum going forward in time under influence of the 
source Jj^ and then returning back in time under the influence of the source J_ (hence the term the closed-time- 
path formalism broadly used in literature [l^). From the path-integral representation of the in-in generating 
functional (for details see e.g. p^). 



one obtains the expectation values upon differentiating with respect to J+ and J_ and then setting 

ip{Xn)}\^,tin), 



(69) 



lSJ+{xi) l6J+{Xn) 
{n, tin\{T(p{yi) ■ ■ ■ (fiiyrn)} X {T<f{xi 



(70) 



with the boundary (initial) conditions for the m-vacuum state, and the constraint for the fields on the future 
boundary ip+{tout) = f-iiout)- In (1701) T and T denote time and antitime ordering operators, i.e., T(p{x)(p{x') = 
e(a;° - x'°)ip{x)(p{x') + e(a;'" - x°)ip{x')ip{x) and Tip{x)ip{x') = e{x'° - x°)ip{x)ip{x') + e{x° - x'°)Lp{x')ip{x). 

The generating functional is now an integral over two fields which we shall conveniently combine into a column 
matrix: 



with the source fields: 



J 



(71) 



To calculate the generating functional Z[J^, J+] is a rather formidable task, which principally can be performed 
only in the perturbative setting. Hence we decompose the action S[(p] into its free (S'i?[(^]) and interacting part 



d x-ipVoif, 



(72) 



where for a nonminimally coupled scalar field the kinetic operator Vq in ((7^ takes the form: Dq = D — — £,R, 
where □ is the d'Alembertian operator, m is the scalar field mass, R is the Ricci scalar and D is the dimension 
of space-time. The in-in generating functional (j69p now becomes 



Z[J-,J+] 

Next we shift the fields 
which results in 

z[J-,J+] -- 



^zS+[-i5/5J+]-iSY[iS/SJ^] ^ 



(73) 



X exp 



d^x 











iJ+ 
-iJ- 



iA_+ zA__ 



-iJ- 



(74) 



(75) 



«Ah 

2A_ 



iA_| 
«A_ 



{x]x') 



{x') 



X exp y d'^xd'^x —{iJ+,—iJ^)(x) 
where the matrix of A's is a 2 x 2 Keldysh matrix of (free) propagators, defined as the inverse of V^s: 

(^)('!a^'' ]^^-^{x,x')^^5^{x-x'). (76) 





-Po 



Here «A_|_+ is the free Feynman (time ordered), zA is the free anti-Feynman (or Dyson) propagator and i/S.^ 

and «A y. are the (positive and negative frequency) free Wightman functions. If we now insert Eq. ()76p into 

Eq. (j70p we obtain relations between expectation values and propagators: 

iA++(a;, x') = {n\T {ip{x)ip{x')) \n) = 8(2;° - a;'")iA_+(x; x') + e(x'° - x")iA+_(a;; a;'), (77) 
iA__(a;,a;') = (17|T (^(a;)(^(x')) = 6(2;° - a;'")iA+_ (x; x') + e(x'° - a;")zA_+(a;; a;') ■ (78) 
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It is worth noting that these relations hold both for free and dressed propagators. When is a massless minimally 
coupled scalar field, i.e. "Dq = d^, the free propagator equations reduce to 

d'^i/^++{x]x') = i5°{x-x'), d'^il::^+-{x;x')^Q , (79) 

92iA__(a;;a;') = -i5^{x-x'), d'^il\^+{x;x') ^ Q . 

The vacuum solution to all four propagators can be written compactly as 

^a..(.;.') = 4J^(^) , (80) 

where Aa;±± are the four Schwinger-Keldysh length functions: 

1^x1+ = \\x-x'f -[\t-t'\-i5f , /^x\_ = \\x~x'\\^ -{t-t' + i5f , (81) 
Aa;^_ = \\x-x'\\^ -{\t-t'\+i5f , 1^x1+ = \\x - x'f - {t - t' - idf , 

where 5 > is an infinitesimal parameter, which defines the time integration contours needed when propagators 
act on functions or distributions. Notice that Ax^^ and Aa;?._ are Lorentz invariant, such that the Feynman 
and anti-Feynman propagators in (|80p are manifestly Lorentz invariant. 



Appendix B: Classical Newtonian potential 

In order to obtain the Einstein-Hilbert action to quadratic order in h^i, (see Sec. II) we first expand the Ricci 
scalar in powers of h^^y. The general result is (see, for example. Appendix A of Ref. |1Q|). 



-gR 



pa ^(7^ 



where 



11 3 

'^{daKa){dphp.i3) - -{dfihpy){dahpa) + {dcrhp,,){dphai:i) + -{dphp^){dyha/3] 



-.9 = 1 + Ih + - hpyh^") + 0{hl,) , .9'^- = T?^"- - h^'^ + h^^hf"^ + 0{hl,) . 



(82) 



(83) 



Notice that (f82|) contains all orders of hp^, and hence can be used to obtain the cubic, quartic and higher order 
gravitational vertices. When Eqs. are combined, we get for the linear and quadratic contribution to the 

Einstein-Hilbert action, 

= (i + ^) [Opdyh^^" - d^h) - hp,{2d''d^h''p - d'^hp'^) + h^-'dpdyh 

- lidphnid'^h^^) - + (9^/i)(a'^ V) + ^id'^h-'nidM + 0{h%) (84) 

= -hkd'^h + ]^hdpdyhP''' ~ ^h'"'dpd„hl + ^h'"'d^hpy + ©(/ij^^) + tot. der. terms , 

where d"^ — i^f^^dpO^, h^"" — rj^^^'r]'"^ hp^, 9^ — rj^^^d^, h — rj^^^hp^. In the last line of ([Ml) we dropped several total 
derivative terms (which contribute as boundary terms, and hence do not contribute to the equations of motion). 
The last line in ((M)) can be used to recover the quadratic action (jH). 

Next, for pedagogical reasons here we present a derivation of the classical potentials ©, which solve Eq. ([S|): 



where — IGttGn and the Lichnerowicz operator is given by 



(85) 



L^^ypa = d(p'q„)(pdy) - ^Vp.{pVa)^d'^ - T^iVpcrdpdj, + i]pi,dpda) + ^-rjp^ijp^d'^ . 



(86) 
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Let us first write the (00), (Oi) and (ij) components of Eq. ([S5 



^{d,d, - S,,V')h,, = yAf<53(f) (87) 



+ \5^A-^l + y^)hu-\{5^J^k^lhkl +^,^Jhu) - 0, (89) 

where — Sijdidj = didi and hij are the spatial components of h^^. 
Next we perform the standard scalar- vector-tensor decomposition of hf^^: 

ho^ = nf + d,a , h,, {d,d, - V^) h + d,hj + d,hj + hj^ , (90) 

such that diuf = = dihf , hj^^ — 0, dihjj^ ~ and Tr[/iij] = h. Consider now an infinitesimal coordinate 
transformation, — > + £,^{x). Then, as a consequence of the tensorial transformation law, the metric 
perturbation on Minkowski background transforms as, 

V(^) KA^) - d^U^) - d^Ux) + Oie,h^,^h) , (91) 

(with — "q^j^vC) where it is also known as the gauge transformation of general relativity. In the light of the 
metric decomposition (|90p . it is also convenient to decompose into two scalars and a transverse vector as 
— {^Ot^I idiC)^ where di^J — 0. With this ((9T|) can be rewritten in components as 

^00 hoo - 2do£,o , nj ^ nf - d^^J , ^ cr - 5oC - , 

hj;" ^ hj;" , hf^hj-^f, h-^h-2C, h^h-2V'C (92) 

From these relations it then immediately follows that the following two scalars, one vector and the tensor are 
gauge invariant: 

hoo - 2doa + dfh , h-V^h, nf - dohj , hj^ . (93) 

This means that out of the ten components of /i^^, six are gauge invariant (and correspond to observables) , 
while the remaining four components are gauge dependent, and have no independent physical meaning. It is now 
instructive to rewrite Eqs. (|87H89| in terms of the scalar, vector and tensor components of /i^i/ defined in (|90l) . 
The resulting equations are 

V^{h-V^h) = -^M5%x) , (94) 

V2(nf - dohf) = , dod^{h - W^h) = , (95) 

id,dj-S,jW^)[^{hoo-2doa + dfh)~^{h~W^h)\ =0, d^{h-V^h)=0, dod,{nJ - dohj) = , (96) 

{~dl+^')hl^ = 0, (97) 

These equations possess a number of remarkable properties. First, Eq. (j97p is the only dynamical equation and 
it is the wave equation for gravitational waves, obeying the standard Lorentz covariant wave operator. Second, 
all equations are gauge invariant (see Eq. ([93])), and hence can be written in the gauge invariant form as 

V^^r = 47rGArA«3(-) ^ = = , (a,aj -5,jV2)(-$ + ^') = , So^^nf = V^nf = O , (98) 

where we introduced the usual Bardeen potentials, 

^ = -\{h~^^h) , $ = -i(/ioo- 290(7 + /i) , (99) 
and a gauge invariant shift vector, 

nf = - d^hj . (100) 
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Eqs. (1^5]) can be easily solved assuming that $ — 5* and nf vanish at spatial infinity: 

GnM 



$ = * = -- 



(101) 



where r — \\x\\ and we made use of the Green function G{x ) — — l/(47rr) for V'^G{x ) — 5^{x ). Notice that there 
are 3 equations for ^ in ((98|). These equations are not redundant. In fact, they assure that the linearized version 
of the BirkhofF theorem holds: the potentials of linearized gravity generated by a static point mass must be time 
independent. Notice finally that, the linearized Eqs. (IMl - f571) do not tell us anything about the gauge variant 
quantities (among which are two scalars and one vector), which is also what one expects. 

Most of the literature does not make use of the (gauge invariant) Bardeen potentials to specify the gravitational 
response to a static point mass, but instead one fixes a gauge. The main purpose of this Appendix, in which we 
show how to construct classical gauge invariant potentials, is to facilitate comparison with the literature, which 
uses a variety of gauges. For example, in the longitudinal (Newton) gauge, the scalar metric sector is fixed to be 



(102) 



In this gauge we have simply 

$ = 



^00 



6 



(no summation) . 



(103) 



Together with poip . these equations imply which is the classical Newtonian potential that we use in this 
work. 

On the other hand. Park and Woodard [iB] make use of the non-diagonal gauge. 



^00 — 



2GnM 



2GnM 



• i' 3 ■ 



(104) 



This gauge can be obtained from the two scalars of hij in (|90|) by observing that for a radial scalar function 
h — h{r), di — fid/dr {fi = x'/r) and didj — {Sij /r){d/dr) + fifj[{d^/dr'^) ~ (l/r)(d/(ir)], from which it follows. 



ij Jscalar 



_(d^ 



1 d 

\ dr^ r dr 



(- 



1 d 
r dr 



Comparing this with (|104p reveals that in that gauge 

h-(— -— 
V dr^ r dr 



2GnM 



To check these relations, note first that, up to a constant, h = —2GNMr. Equation ((99|) then implies 



2 

\dr-^ 



2 d 
r dr 



GnM 



(105) 



(106) 



(107) 



which agrees with ()10ip . implying that Eqs. ()104|) represent a correct classical solution. 



Appendix C: Expanding the bubble diagram 

In this Appendix we expand the scalar field bubble diagram in Eq. ()17p shown in Fig. [3] in powers of small 
(quantum) metric perturbations Sg^"'{x) = -/i^''(a;) + h^'^'' {x)h^,{x) + 0{{h^"'f). 
Let us begin by noting that the scalar propagator in (1171) obeys the equation 

{-9^'\x)f/^U'i^I^'^\x,x', [9%,]) = 9,g(^)"^(x)[-.gW(x)]V2a,A-''(x;x'; = {a^yH^ix - x' ) , (108) 



such that y/-g('>){x)ai '6'^{x - y){(J^Y^ is the operator inverse of A''''(2/; x'] [g^^)p,]). <J^ = diag(l, -1) in pllg)) is 
the Pauli matrix. 
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Writing gjj'} — + 5(7^^, we can expand the bubble diagram in (ITTl) in powers of 5gf^^. To proceed, let us vary 
the first bubble diagram in ([T7)) with respect to g^''{x) (for notational simplicity we drop the Keldysh indices 
a, 6, .. = ± and the background index (&)): 



- -TvH/^p[y-y-[g%])] 



(109) 



d^y / d'^zdy 



where Ap{z;y; [ga'p']) = A++(z;i/; [ga'p']) denotes the Feynman propagator. Now, because of the 5^(y — x) 
function, all of the ?/'s in the square brackets can be replaced by x's, and a part can be pulled out of the integral 
to obtain 

{-\V^)9''^{^)9^''{^) + V^)SiX))l f d^'y f d^zS''{y~x)[d^pS^{x~z)]Ap{z;y; [.g^,^) [a^<5^(y - x)] 



Next, the derivatives 8^6^ {x — z) = —dp6^{x — z) and 9^ can be moved by partial integration to act on the 
propagator Af(z; y[gap\), and finally the two integrals can be performed on the expense of the J-functions. The 
result is 



5g^^''{x) 



-Trln[AF(y;2/; [ga-p'])] 



{\^J^9iF)9''^{^)9,u{x)-^/^)5l5l^) [dldliAF{x-y- [g^,^,])] JUO) 



1/1 
2 



This is the general expression for the tadpole diagram, i.e. the linear coupling of the stress energy tensor to 
metric perturbations. 

In order to get the vacuum polarization we need to vary (jllOp one more time with respect to g'"'{x'). To 
facilitate the variation procedure, recall the operator inverse identity: 



rf^z< V^^5^^(z')9|'<5^(z' - z)AF{z-y- [g^,p,]) = 5° {z' - y) . 
Varying this with respect to g''"'{x') yields 



d^z / d^z 



'Af{x;z'; [<?„'^'])< V^^.9^^(z')3p.(^')'5^(^' - + ^^^)SlSi^S^ {z' - x'))dl' (z' - z) 



X AFiz;y;[ga'p']) + / d^zd^{x-z) 



S 



—AF{z;y; [ga'fi']) = 0. 



SgP<^{x' 

Setting y ^ X and performing similar steps as above, this expression can be recast as 
—AF{x;y; [ga'p']) 



^{x') 



\,/^^)g'^\x')g,^{x') ^r^)Sl6i^) 



d^'Apix; x'\ [ga'/3']) Af(x'\ x\ [ffa'^j']) 



(111) 



Inserting this into (|110p we get for the second variation of the bubble diagram: 



bgP^ix') Sgf^^ix) 



-Trln[AF(2/;y; [ga'fs'])] 



^y=^(ig"'^(x)g,.(2:) - S^X)) [d:d^,'^AF{x;x'; [g^,^,]) 



dpi\AF{x';x; [g^,^,])] ^/^:^)(lg-f\x')g,,{x') - SJ/.^) 



-rl 



\/-9ix) -^9paix)g"l^{x)gf,^{x) - -SfpS^^g^,^{x) - -Sf^^gp^ix) + -g"'^(a;)g^(p(x)g^)^(x) 



X S°{x ~ x') d^d^p iAf{x; x'; [ga'p']) 



(112) 



The last two lines come from varying the metric tensors in (jllOp . Putting everything together, we get the following 
expansion of the scalar bubble diagram: 



lTrln[A^^(y;y;[5S]) 



-Trln[AF(2/;y; [ga'p'])] 



d''x^/^^Sg^-ix)Tp,ix; [g^,p,]) 



(113) 



i / d^xV^^) [ d^x'V^^(5g^''(x)[^,n,.](x;x';[g„,;3'])<5ff'''^(a;')+0(('55"^)'), 
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where 



2 V2 



X' 

5 



(114) 



d^d^'iApix; x'; [5a'/3']) OpfiApix'; x; [^a'/J']) 



1 



j9pT{x)g°'^{x)g,,u{x) - ^6"pS^^gp^{x) 
d^dpAFix;x';[gc,'0']) . 



^{,.^^)9p'^ix) + j9'"^ix)gt,ip{x)ga)uix) 



(5^(a;-x') 



\/-9{x) 



(115) 



The interpretation of this result is simple. The first term in (|113l) is the scalar IPI bubble in a classical background 
9ap, the second term is the IPI tadpole, with the metric perturbation Sg^'^{x) — -~h'^'^{x) = —g'^°'{x)g'^^ {x)hap{x) 
and the matter stress energy tensor T^^, given in (|114p . The last term in (|113p is the contribution that is 
quadratic in the metric perturbation 5g^'^ ~ —h^^^ with ^u^pa{x',x' ;[gai p']) being (a part of) the graviton 
vacuum polarization (calculated on a gap{x) background), whose explicit form is given in (jllSp . The first two 
lines in (jll5p represent the nonlocal part of the graviton vacuum polarization shown by the middle diagram 
in Fig. [21 and they perfectly agree with the graviton vacuum polarization shown in the first line of Eq. (I119p 
in Sec. IIIII (■f — J> and gap — >■ rjap)- The last two lines in (|115l) represent a part of the local contribution to 
the vacuum polarization. The other part is obtained by observing that the quadratic term in the expansion 
Sg^^ix) = -h'"'{x) + h^'p{x)hP'' {x) + 0{{h^''' {x)f) in the tadpole in (fm]) yields the local contribution to the 
vacuum polarization of the form 

5[p.n,,]{x-x'-[g^.p>]) - [-]^g'^^{x)g,(pg,)p{x) + 5'^^g,)f^/^^]5^ . (116) 

such that 

'][9a'P']) = [p^Ii-pcr\{x\x';[ga'l3'])+5[^,^Iipcr\{x\x'\[ga'p']) (117) 

represents the complete graviton vacuum polarization, which correctly includes also the local contribution from 
the quartic interaction (1131) (in the limit when g^^ — > r/^^ and ^ — > 0). A technique that is equivalent to expanding 
the 2PI bubble diagram around the Minkowski background was used to calculate the leading quantum correction 
to the Newtonian potential of a static point particle in Refs. [26, 2^. In the limit when time transients are gone, 
our results presented in Sec. IIVI indeed agree with those of Refs. [1^ [13], just as one would expect from the 
analysis presented in this Appendix. 

And finally one more remark. From Eqs. (jll3Hll4)) and the preceding analysis it follows that the nonlocal part 
of the graviton vacuum polarization is obtained by 



{[pu^pa] (x;x')}no„_ 



local 



\j9{x') Sgpa{x') 



[Tp„ix; [ga'P'])] , 



(118) 



where the tilde on the functional derivative (d) signifies that it acts on the propagator(s) in T^^ only. A special 
case of this relation can be found in Eqs. (12-15) and (47) of Ref. [lH, where the non-local part of the graviton 
vacuum polarization tensor around the Minkowski vacuum is expressed in terms of the non-local. Wick contracted, 
stress-energy-stress-energy correlator . 



Appendix D: Non-local contributions 



In this Appendix we carry out thorough calculation of the graviton vacuum polarization. By Wick-contracting 
all non-coincident scalars and making use of the propagator definitions (|77H78p , the graviton vacuum polarization 
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inferred from 



becomes 



1 



4C 



- v'^'^'vpcrdad[yiA±±{x; x')dg,^dpiA±±{x; x') + 9(a^A±±(x; x')dfj)D'p^iA±±{x; x') 
d'a'd(p^iA±±{x;x')d^)d'p,iA±±{x;x') - Tf^-qp,^d'^,d(^^iA±±{x]x')ds)d'p,iA±±{x]x') 



+ d[^,tA±±{x;x')d'^,^D^,tA±±{x;x') ( ^v^^'^'vpa " 4"'-^? 



9^9(^«A±±(a::;a;')c}^)9^zA±±(a;;a;') - rj'^ * 77p^a^a(^,iA±±(x; x')9^,)9t.zA±±(a:; x') 

+ 9(^«A±±(a;; a;')9^)D^^7A±±(a;; a;') - rf^'q^^d(^^iA±±{x\ x')ds)D'p^iA±±{x; x') 
+ d[piA±±{x; x')d'„-)Df,,yiA±±{x; x') - rf' ^ ripad[^,iA±±{x;x)d's,)D^uiA±±{x]x') 

+ \i^±±{x\ x')Dp^D'p^iA±±{x; x') + ^D^^iA±±{x; x')D'p„iA±±{x] x') 



(119) 



The first line in (jll9l) corresponds to the graviton vacuum polarization tensor of a minimally coupled massless 
scalar field already calculated in ^16]. Our analysis in Appendix C shows that this expression can be also 
obtained by expanding the 2PI bubble diagram around the Minkowski space, i.e. take Eq. (jllSp in the limit 
when Qfj^i, — r/py. The vacuum polarization (|119p contains a local divergence. In order to extract it, the following 
identities can be employed ^31], which are easily derived from (150115^ : 



d'pd^iA++{x-x') = 5l5li5''{x-x') + 
d'pdaiA--{x]x') ~ 



n , r(f) 



d'pdaiA^+{x]x') 
d'odaiA+^{x;x') 



5l6l^5^{x-x')^^^ 



D 



AXaAxi3 



(Ax^+)^/2 (Ax^+)(^/2)+i 

AXaAXjS 



- D- 



27r^/2 



D- 



_{Axl_)D/^ {Axl_)(D/2)+i 
AxaAxp 



(Ax^+)^/2 (Ax^+)(^/2)+i 



D- 



Axq Ax^ 



(Ax2__)^/2 (Ax2__)(-D/2) + l 



(120) 



In order to reduce the i^-dependent terms we also need the following derivatives: 
^,^'p,^',,^A±±ix; x') = (<50<5° 9^, + 6151,8'^, - S^, S^, 8^)16'' {x - x') 



D 



r 



27r^/2 



rip^a'Ax^ + -qp^aAXa' + 1 ]^'aAXp, _ ^jj ^^AXp'AXa'AXa 



Ax 



D+2 

±± 



dpd.d'p,d',,iA±±{x;x') = {5X'd'.'dp+S'.,S°d.d'p,+S°X,d'^,dp + 5^p,S°^d.d',,-5"^^^^^^^ 



D 



r(f) 

27r^/2 



AXpAXa-AXp' AXa-' 



Vp'a'Vpa+Vp'aVcr'p + ria'aVp'p , , ^^^^ , 



{D + 2) 



±± ^■^±± 

rjp'a'AXpAXa- + r]p>crAXpAx^' + r]pipAXa'AXa + riper' Ax p'AXa + -qp^AXpi AXa> + r]aa'AXpAXp, 



Ax 



D+A 

±± 



.(121) 



When these expressions are inserted into the graviton vacuum polarization (|119p . the terms in (jl20H12l|) that 
contain delta functions do not contribute in dimensional regularization, since they hit a propagator which contains 
a function of the type (Aa:)~", where a — a{D) is a D-dependent power. One can show that rigorously by the 
analogous reasoning as was used in Sec. IIII Al 

The next step in the reduction of the graviton vacuum polarization tensor is the extraction of derivatives. In 
the following we express various terms appearing in the vacuum polarization as an operator acting on the single 
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function, 1/Ax (which stih contains the local divergence that we are after): 







1 




[AD{D-l){D-2)^ } 




Axf^Axi, 






Axll+' ~ 1 


[81)2 (£»-!)(£)- 2)2 


} Axil-' 


Ax Ax ^ Ax p Ax a- 







At2^+4 



\QD^{D + 1){D -!)(£>- 2)2 ' 
^ 1QD{D + 1){D -1){D ~2Y 



(122) 
(123) 
(124) 



At2^- 



Before reducing all the terms of the tensor (|119p it is instructive to note the structure of the parts linear and 
quadratic in ^ which are basically the same up to a multiplicative constant. To be more precise, in the second 
and third terms that are linear in ^ the following operator can be extracted: 



-{D^-D- 2)D^,Dp,, 
while in the term that is quadratic in ^ we can extract 



4(^2 - l)Dp,D^ 



pa- 



(125) 
(126) 



After some algebra, and after combining the linear and quadratic contributions in ^, the graviton vacuum polar- 
ization (I119P can be written as 



z[±,n±] ix;x') 



±K 

~2" 



±K 

~2~ 



r2(f) 



16nD{D^ - l)iD-2) 



(127) 



Upon performing one more partial integration: 

1 



D 



D- 1 



D 



fii^ D pa 



Ax 



2D-4 
±± 



92 



2D-4 



2iD-S){D-4) A: 



2D-6 



(128) 



we see that {D — 4) appears in the denominator, signifying a divergence. This divergence can be removed by 
recalling the identities (that were employed when constructing the propagators ((80|) ). 



1 



Ax 



D-2 

++ 

1 



4^^/2 

r(f-i; 



iS"{x-x'), 



Ax^-2 



Based on these identities and Eq. (jl28l) we easily obtain 



{x — x') , 



92 



Ax?-' 



1 



D-2 



0, 

: 0. 



(129) 



1 



52 



1 



D-4 



A; 



2D-4 



2(D-3)(D-4) \Ax 



2D-6 
±± 



A: 



D-2 
'±± . 



,3\±± 



2^D/2^I3-4 



r{§-l)iD-3){D-4) 



iS^{x-x'), (130) 



where a — diag(l,— 1) is the Pauli matrix and fi is the mass scale appearing for dimensional reasons which 
signifies the renormalization scale. The identity p30p is judiciously constructed, such that (Axj__|_)2^-^ has been 
split into a piece that is finite in Z? = 4 and a local divergent piece. Indeed, when expanded around D = 4, 
Eq. (|130p can be recast as 



2D-8g2 



Ax 



2D-4 
±± 



32 



[ln2(^2^^2^^) _ 2 \n{^i^Axl±) + 0{D - 4)] 

2^^/2/i^-4 



r(f ~l)(i?-3)(i?-4) 



i5"{x-x'), 



(131) 
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where we made use of 



,15-4 



Ax 



±± 



,2_D-8 r 



D 



■In (Ai'A4±) + 0((i^-4)2) 



and 



ln(^^Ax^±) ^ 1^2 



ln2(^2^x^^) - 21n(^"A.T^±) 



(132) 



(133) 



With the resuhs outlined in this Appendix it is possible now to express a graviton polarization tensor in terms 
of the sum of the non-local finite and local divergent parts ([27]). 
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